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Abstract 

We calculate all the experimentally relevant branching ratios of the Higgs 
bosons of the Minimal Supersymmetric Standard Model, paying particular 
attention to the contributions from below-threshold decays. We show that 
in some cases these can significantly change the pattern of branching ratios 
calculated without taking off-shell effects into account. 
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1 Introduction 



It is likely that the first experimental evidence for supersymmetry will come from studies 
of the Higgs sector. In supersymmetric extensions of the Standard Model (SAi) the Higgs 
sector invariably has a rich structure, which in the case of the Minimal Supersymmetric 
Standard Model (AiSSJ^) is also highly constrained (see for example Ref. []I|). Finding 
one or more of these Higgs bosons is one of the major goals of present and future high- 
energy colliders. 

Many years of phenomenological studies P, H, ^ || have produced several 

important lessons for Higgs searches. The first is that almost all decay channels have 
large SAi backgrounds. It is therefore important to have very precise calculations for 
the Higgs production cross sections and branching ratios, so that small signals can be 
unambiguously identified. A second lesson is that 'below-threshold' decays can be very 
important. By this we mean Higgs decays via intermediate states in which one or more 
particles are off-mass-shell. A simple example of this is provided by the SA4 Higgs decay 
—* Z^Z^ — i> 4/^. In the 'on-shelF approach, one would calculate the branching ratio for 
this as BR{(f) Z^Z^) x BR{Z° l^l'?, with ^ Z^Z^) = for < 2M^o. 

This, however, is too naive. The below-threshold 'off-shell' decay (p Z^{Z^)* —>■ 41"^ is 
non-negligible, and in fact provides a very important Higgs signature for 0(130 GeV) < 
< 2Mzo. 

Although below-threshold decays have been extensively studied (see below) for the 
SJ^, to date there has been no equivalent study of the J^SSAi. In fact all J^SSAi Higgs 
phenomenological studies have so far used branching ratios based on on-shell calculations 
(except for the trivial cases of H ^ (W^)* , Z'^ (Z^)* which are easily obtained from 
the corresponding SAi results). In this paper we present a complete calculation]^ of all 
relevant AiSSAi Higgs branching ratios, paying particular attention to below-threshold 
contributions. Of course it is impossible to cover all regions of the AiSSAt parameter 
space. Our aim is to identify those situations where below-threshold decays are important, 
and those where simple on-shell calculations are sufficient. We illustrate our results by 
numerical calculations using 'typical' parameter values. Our aim is to provide the necessary 
theoretical tools for precision AiSSAi Higgs phenomenology at future colliders. 

The paper is organized as follows. In the next section we give a comprehensive review 

^For simplicity we ignore all direct decays involving supersymmetric particles 
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of SAi and AiSSAi branching ratios. In Section 3 we present numerical results for 
some typical A4SSA4 scenarios, highlighting the differences between off-shell and on- 
shell calculations. In Section 4 we present our conclusions. 



2 SA4 and A4SSA4 Higgs Branching Ratios 



The branching ratios of the A4SSAi Higgs bosons have been studied in various papers. 
A good review of the early works on this subject can be found in Ref. where all the 
most relevant formulae for on-shell decays have been summarized. These include the 
iS^Vl-like two-body decay^ $ — > // (where / represents a generic massive fermion) 0, 
H — > W^W^ , Z^Z'^ |T0| at tree-level, and the one-loop induced decays $ Z°7,77 
TT| and gg fl^, generalized to the MSSM [l^, IT 



and the 55 -specific decays 



H 



W^K H hKAA, $ 



xfxj (chargino decays) and 



^ ~^ XiXj (neutralino decays) [0. At the time of Ref. |jT| only the tree-level Higgs mass 
relations had been computed, and the on-shell decays h AA and H H^H^, Z^A 
were kinematically forbidden. We now know that they are allowed at one-loop and the 
corresponding formulae can be easily obtained, for example from Ref. Finally, the 

one-loop induced decays of the charged Higgses H"^ W^Z^, W'^'j have also been studied 

Recently, higher-order corrections to most of these processes have been computed. For 



the SAi, one can find in the literature the QCD [|l6l and EW [p^Tj p!8| , |19|, ^ radiative 
corrections, and their interplay ||21|, ^ for the //-channels. The EW corrections to the 
W^W^, Z^Z^ decay rates 1T|, g|, |4[, and both the QCD |2|, ^ and EW corrections 



to the 2'°7, 77, gg decay rates are also now available. A detailed and updated review of 
the higher-order corrections computed to date within the SM can be found in Ref. 



As at tree-level, the higher-order SM. corrections can easily be extended to the 
MSSM case. For the A^55Ai-specific processes, both the QCD |2|| and the EW p9 
corrections to the fermionic decays — > //' of the charged Higgs bosons, and the EW 
corrections |Q to the decays involving neutral bosons (such as H ^ hh and A Z^h) 
have been computed. 

Within the MSSM, a systematic analysis of all the interesting decays of both neutral 



^Throughout this study we use $ to denote the set of neutral AiSSAi Higgs bosons h, H and A 
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and charged Higgs bosons, with reference to the search strategies at the LHC proton- 
proton colhder and taking certain of the above results into account, has been presented 
in Ref. using tree-level formulae for masses and couplings. This analysis has been 



updated in Ref. [31 1 to include new results for the one-loop corrections to the latter |32, 33 



Similar studies have been reported in Refs. and extended to the LEP I and 

LEP II e^e~ colliders in Ref. [^]. For a NLC e^e~ collider the corresponding analysis has 
been performed in Ref. PS|. 



Apart from the processes H —> W W^*, Z Z [39|, the results on decay widths and 



branching ratios reported in Refs. [^T], ^ ^ ^ ^ ^ have been obtained for the case 
of above-threshold decays only. It is our aim in this paper to generalize these results to 
include also below-threshold decays. In particular, we study the two-body decay channels 



^ ^ SS,CC,bb,t*t* , fJ.'^^ ,T^T , 

$ ^ ■^±*yf/T*^ Z°*7, 77, gg, 

H h*h\ H A*A\ A Z^*h\ 
for the neutral Higgses ^ = H, h and A, and 

—y cs, t*b, fi^u^, t^Vt-i 
for the charged Higgses if^'s. We do not consider (i) the tree-level decays 



H Z^*A* 



h A*A\ 



(1) 

(2) 
(3) 

(4) 
(5) 

(6) 



since these are only possible in regions of the AiSSAi (M^, tan (3) parameter space already 
excluded by LEP I (for nit < 200 GeV) pOI, and (ii) the one-loop induced decays 



iy+*7, 



(7) 



since these have very small branching ratios |15 



In calculating the decay rates for the processes (|1|)-(|^) we have adopted the off-shell 
decay formulae presented in Ref. for $ 1^ W*, and in Refs. g3[ for $ ^ 

Z°*7. For the other cases not involving decays to top quarks we have recalculated the off- 
shell formulae at tree-level, tracing the corresponding amplitudes squared and integrating 



3 



analytically over the appropriate off-shell phase space. We then integrate numerically 
over all the possible virtualities of the final-state particles for all the processes. This is in 
contrast to the procedure sometimes used of constraining one of the decay products to be 
on-shell when the decaying Higgs boson mass exceeds its rest mass. Although this allows 
the width to be computed analytically, it can give misleading results near threshold. Our 
procedure is only to use such 'integrated' analytic expressions well above threshold. We 
therefore proceed in the following way. With A the decaying particle and B^*^ and C^*^ 
the decay products, which can be either on-shell {B, C) or off-shell (i?*, C*), we define as 
the decay width of the channel A the function 

^ ^ J Jr AiQl-Mlf + iaBTsf] 



2 

xr{A^B*C* 



dQ'cMcTc 



where Mb{c) and '^b{c) are the mass and the width of the particle B{C) with four- 
momentum squared Q%{c)i respectively, and T{A — > B*C*) the off-mass-shell decay width 
into B and C. The two-dimensional integration region R is defined by 

Qb > nib, Qc > m,, Qb + Qc< Ma, (9) 

where mb{c) represents the sum of the rest masses of the decay products of the particle 
B{C). The change of variables K 



,2 (g|-Mi)2 + a|.^ 



Qj,~M'^=MxTxtanex, ^ dQj, = ^ dOx, X = B,C, (10) 

then gives an integrand smoothly dependent on the integration variables. For the numerical 
evaluation we use VEGAS |45|. 



The expression given in (H) is a good approximation when Mx ^ fn^ {X = B,C) 
so that Tx oc Mx- We therefore use it for decays to vector- vector, vector-scalar and 
scalar-scalar pairs. It is not however a good approximation for decays involving the top 
quark, e.g. $ ^ or tb. For such decays we must perform an exact matrix element 

calculation for the complete process, e.g. $ — > bW~^bW~ or —>■ bW~^b. 

We are not interested here in studying in detail the dependence of the processes (|l])-(|I) 
on all possible higher-order corrections to the widths and branching ratios. Instead we 
wish to find out which below-threshold AiSSM. Higgs decays are important. We therefore 
include only the larger non-SUSY QCD corrections (when known), and not the EW ones. 
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In any case, our procedure can readily be adapted to include such corrections if necessary. 
We have included the modifications due to QCD loops in the computations following 
Ref. |]l6l for the processes $ qq, Ref. |^ for if qq', Ref. ||2^ for h — > 77, Z°7 



and, finally, Ref. ||26| for h gg. In particular, since the QCD corrections to the top 
loops in ^^7, 77 and gg decays, as given in Refs. |2^, are valid for M$ < 2mt, we 
have implemented these only for the lightest AiSSAi neutral Higgs, i.e. ^ = h. We do 
not include any threshold effects due to the possible formation of tt-bound states in the 



one-loop induced processes, when the Higgs mass M<j, reaches a value ~ 2mt |46|. 

There is a slight subtlety concerning the 'running' of the fermion masses. For the 
decays involving light quarks $ — gg (g = s,c,b), the use of the running quark mass 
^qiQ"^ = ^1) takes into account large logarithmic corrections at higher orders in QCD 
perturbation theory [^]. In principle, one could imagine using the same procedure for 
$ — > in the limit M$ ^ m^. In practice, however, we are interested in the case of 
Mis,/mt ~ For the top-loop mediated decay $ gg, it is known that the higher- 

order QCD corrections are minimized if the quark mass is defined at the pole of the 
propagator, i.e. m'^{Q^ = ml) ||4^. For the pseudoscalar^liggs, the decay rate for gg 
has a local maximum at M4 = 2mt. To be consistent, therefore, we must use the same 
top mass {mf{Q'^ = mf)) in the decay rate for $ tt. 

To avoid a proliferation of additional parameters, we neglect small chargino and neu- 
tralino contributions in the $ ^°7, 77 decays. To further simplify the discussion, we 
assume a universal soft supersymmetry-breaking mass ||32|, ^ 



2222 /1 1 \ 

mQ=mU = ^D=^q^ (11) 

and negligible mixing in the stop and sbottom mass matrices, 

At = Ab = fi = 0. (12) 

Moreover, if we also neglect the 6-quark mass in the formulae of Refs. ||3^, the one- 
loop corrected masses of the AiSSAi neutral CP-even Higgs bosons can be expressed in 
terms of a single parameter e given by 

In f 1 + (13) 



STT^M^i sin^ Ow 



mi 



Diagonalization of the mass-squared matrix leads to the expressions 

MIh = ^[Ml + M|o + e/sin2/5] 
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± {[{Ml - M|o) COS2/3 + e/sm2/5]2 + (M^ + M|o)W2/3}'^' , (14) 

while the mixing angle a in the CP-even sector is defined at one-loop by 

^ ^ (Ml + M|„)sin2/5 ^ n .^ 

tan 2a = ^ ^ ^„ ^ ^ — — < a < 0. (15) 

(M|-M|o)cos2/3 + e/sin2/5' 2 " ^ ^ 

For the AiSSAi charged Higgs masses we have maintained the tree-level relations 

M|±=Ml + M^±, (16) 

since one-loop corrections are quite small in comparison to those for the neutral Higgs 
bosons PB| . 

For the above choice of parameters we also use the leading one-loop corrected expres- 
sions for the Hhh and HAA trilinear couplings, which enter in the Higgs branching ratios 
studied here, i.e. 

^Hhh = ^%hh + ^^Hhh, (17) 
^HAA = ^^HAA + '^^HAA, (18) 

with the tree-level relations 

\%fih — ~7r' — sin(/5 + a) sin 2a — cos(/5 + a) cos 2a;], (19) 

2 cos uy/ 

>^HAA = cos2/3cos(/3 + a), (20) 

Z COS u\y 

and the one-loop corrections 



igMzo 3g'^ cos^ 6w cos^ a sin " "^t [^3 ^ m"^ + ^ m? 
2 cos 9w Stt^ sin^ P f^w^ \ ^ "^i + 



i^M^o cos2 sin a cos^ /3 m? + 

A A;^ = - 7^ n ^2 — ^ — ' (22) 

2 cos c^vi/ oTT"^ sm p ^^4/± 



as given in Ref. |31 



In the numerical calculations presented in the following section we adopt the following 
values for the electromagnetic coupling constant and the weak mixing angle, agm = 1/ 128 
and sin^ 9w = 0.23, respectively. The strong coupling constant as, which appears at 
leading order in the gg decay widths and also enters via the QCD corrections, has been 
evaluated at two loops, with = 190 MeV, and with the number of active fiavours Nf 
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(and the corresponding Aj^'', calculated according to the prescription of Ref. |^) and 
scale chosen appropriately for the decay in question, i.e. adopting the prescriptions of 
Refs. m, H. 

For the gauge boson masses and widths we take M^o = 91.1 GeV, F^o = 2.5 GeV, 
Mw± = Mzo cos 6w and Tw± = 2.2 GeV, while for the fermion masses we take = 0.105 
GeV, nir = 1.78 GeV, m, = 0.3 GeV, = 1.4 GeV, = 4.25 GeV and = 175 GeV 



with all widths equal to zero except for F^. We calculate this at tree-level within the 
M.SSJ^, using the expressions (for the above values of and m^) fST 



F(t ^ bH^ 
T{t bW- 



A(M2,„m2,m?)V2 
[ml + ml — M^±)(m^cotan^/3 + m^tan^/5) + Am 



M^±{ml + ml 



m 



2\2 



(23) 



and 51 



F(t ^ bW~ 



|\/,,P^A(M^„m^,m?)V2x 



8V27t' 



l-(^ 

rrit 



l+(^ 

nit 



\ mt J 



\ mt J 



(24) 



where Vtb is the Cabibbo-Kobayashi-Maskawa mixing angle (here set equal to 1), Gp = 
\^g'^ /8M^± is the electroweak Fermi constant, with g = e/ sin^vy and — e the electron 
charge, and A^^^ is the usual kinematic factor 



A(M„ Mb, M,)i/2 _ ^ ^2 ^ ^2 _ 2^^jvf, - 2M,Me - 2MbM,Y 



/2 



(25) 



The first generation of fermions and all neutrinos are taken to be massless, i.e. m„ = 
md = me = m^^ = and m,y^^ = m^^ = 0, with zero decay widths. Finally, the universal 
supersymmetry-breaking squark mass is in the numerical analysis to be mg = 1 TeV. 



3 Numerical results 

It is impractical to cover all possible regions of the AiSSM. parameter space. We choose a 
representative set of figures which (a) describes the decay channels and mass ranges which 
are relevant to experiment, (b) illustrates the importance of below-threshold decays, and 
(c) allows a comparison with previous studies. Thus Figs. 1-4 show the branching ratios 
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for the h, H, A and AiSSAi Higgs bosons respectively, as a function of the mass of the 
decaying particle. Following Ref. we choose two representative values, tan/3 = 1.5, 30, 
which give significant differences in certain channels. In comparing our figures with those 



of Ref. ||3T|, it should be remembered that we have chosen a top quark mass (m^ = 175 
GeV) consistent with the recent CDF measurement. 

Figure 1 shows the branching ratios of the lightest neutral scalar AiSSAi Higgs boson 
h a.s a. function of M^. Note that for this value of m^, the maximum values of Mh are 
98(129) GeV for tan (3 = 1.5(30), achieved in the limit —>■ oo. As expected, the bb and 
r+r^ channels dominate over essentially the complete mass range. The cc and ti branching 
ratios are smaller for larger tan (3, reflecting the dependence of the hff couplings on the 
angles (a, /?). 



The new feature of Fig. 1 compared with the corresponding figure (Fig. 10) of Ref. |^ 



is the appearance of the below-threshold decays h — > W~^W~ , Z'^Z^, especially for smaller 
tan (3 values. In fact for tan (3 = 1.5, the W~^W~ branching ratio is larger than that for 77 
over a sizeable portion of the h mass range (when ~ 88 GeV). Fortunately, because of 
their negligible contribution to the total width compared to the dominant modes, the two 
new channels do not depress the 77 branching ratio, which is very important for Higgs 
searches at the LHC proton-proton collider. This does happen at tan/5 = 30 (and also 
affects the bb and t~^t~ channels), but only in a very narrow mass window, and so is 
unlikely to be phenomelogically relevant. Note that the quasi-singular behaviour of the 
branching ratios near Mh = M™^^ is caused by the mapping of a large range of M4 masses 
into a small range of Mh masses, see Ref. ||3l|] for a fuller discussion. 

In Fig. 2 we show the corresponding set of branching ratios for the heavier of the two 
neutral scalar AiSSAi Higgs bosons H. This has a rather complicated structure which 
is discussed in Ref. [^. In summary, for large tan/? the only important decays are to bb 
and r+r", apart from a very small mass region around Mh = M]^™ where other decay 
channels are important. For smaller tan/3, either the W~^W~, hh or tt channels dominate, 
depending on Mh- The inclusion of below-threshold decays has little effect here, serving 
mainly to generate a smooth transition around the ti threshold region. 

Figure 3 shows the branching ratios for the pseudoscalar A Higgs boson. Here the 
difference between large and small tan (3 is again very marked. For the former, the bb and 
r+r^ decays are completely dominant, and below-threshold decay effects are negligible. 
For the latter, the ti decay dominates above threshold and the Z^h decay channel can 
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also be important. The point which we wish to stress is that the below-threshold ti decay 
plays an important role. To see this more clearly, Fig. 4 shows the corresponding set of 
branching ratios when below-threshold decays are forbidden. (This is the analogue of 
Fig. 12 in Ref. |31|.) In this case for tan/5 = 1.5 there are three distinct mass regions: 
Ma < (Mz + Mh), {Mz + M^) < Ma < 2mt and 2mt < Ma where the dominant 
decay modes are bb, Z^h and tt respectively. However, the middle mass region is strongly 
affected by the below-threshold ti decay. In particular, for 300 GeV < Ma < 350 GeV 
the bb and Z^h decay modes are suppressed far below their 'naive' values shown in Fig. 4. 
This clearly has important phenomenological implications, especially since the leading 
A ^ ti ^ bW^bW~ channel is not readily observable. 

Another interesting effect is seen in the charged Higgs branching ratios. Fig. 5. At 
large tan/5 the only relevant decays are tb and tv. Here the inclusion of the below- 
threshold decay for the former simply gives a smoother transition between the two; compare 
Fig. 14(b) of Ref. At smaller tan (3 there is a window (assuming the CDF central value 
for TJit) for the Wh channel to be important. In fact this would be the dominant channel 
over a sizeable range of Ma in the absence of the below-threshold tb decay. However, 
we see that when the threshold behaviours are correctly taken into account the tb decay 
rate is almost always larger than that for the the Wh channel. Again, this has important 
implications for charged Higgs searches at future colliders. 

4 Conclusions 

In summary, we have performed a comprehensive study of the phenomenologically relevant 
decay rates of the J^SSM. Higgs bosons, paying particular attention to below-threshold 
decay channels. We have discovered several cases where these may be important. For 
tan (3 not too large, there are non-neglible branching ratios for the decay of the lightest 
neutral Higgs boson h into off-shell W^W~ and Z^Z^ . At the upper end of the M^ mass 
range these can reach 0{1Q~^ — 10^^) and (in the case of can be larger than the 

important 77 branching ratio. 

For the decays of heavier neutral bosons H and A, the main impact is from below- 
threshold ti decays. Again for modest values of tan (3 these can suppress the branching 
ratios of phenomenologically important channels such as H ^ hh,W~^W~ , Z^Z'^,bb and 
A Zh, bb, T^T~ below their values calculated without taking the off-shell effects into 
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account (compare Figs. 3 and 4). In the same way, the below-threshold decay — > tb of 
the charged Higgs boson can suppress the branching ratio of the important — > W^h 
channel. 

As the physics studies for future high-energy coUiders become more focused, it is im- 
portant that all the ingredients of 'new physics' search strategies are calculated with as 
high a precision as possible. We beheve that our study of M.SSM. Higgs boson decays 
provides a significant improvement in this respect. 
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Figure captions 

[1] Branching ratios for the lightest neutral scalar AiSSAi Higgs boson /i as a function 
of Mh, for tan/3 = 1.5 and 30. Other parameter values are given in the text. 
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[2] Branching ratios for the heavier of the two neutral scalar M.SSM. Higgs bosons H 
as a function of Mh, for tan/? = 1.5 and 30. Other parameter values are given in 
the text. 



[3] Branching ratios for the pseudoscalar J\4SSA4 Higgs bosons A as a function of Ma, 
for tan/3 = 1.5 and 30. Other parameter values are given in the text. 

[4] As for Fig. 3, but with the below-threshold decays omitted. 

[5] Branching ratios for the charged M.SS A4 Higgs bosons as a function of Mh±, 
for tan/3 = 1.5 and 30. Other parameter values are given in the text. 
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